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Abstract 

We consider two models of one-dimensional discrete random Schrodinger operators 

"00 = V'n+i = in the cases Vk = auok/V^ and Vk = aujk/yk- Here ujk are independent random 
variables with mean and variance 1. 

We show that the eigenvectors are delocalized and the transfer matrix evolution has a 
scaling limit given by a stochastic differential equation. In both cases, eigenvalues near a fixed 
bulk energy E have a point process limit. We give bounds on the eigenvalue repulsion, large 
gap probability, identify the limiting intensity and provide a central limit theorem. 

In the second model, the limiting processes are the same as the point processes obtained 
as the bulk scaling limits of the /3-ensembles of random matrix theory. In the first model, the 
eigenvalue repulsion is much stronger. 

1 Introduction 

We consider two models of one-dimensional discrete random Schrodinger operators given by the 

matrix 
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in the following two cases, referred to as the critical model and decaying model respectively: 

Vk = auk/y/n, Vk = aUk/Vk. (2) 



Here Uk are independent random variables with mean 0, variance 1 and bounded third absolute 
moment. 

We show that the eigenvectors are delocalized and the transfer matrix evolution has a scaling 
limit given by a stochastic differential equation. We show that in both cases eigenvalues near a 
fixed bulk energy E have a point process limit. 

We analyze the limiting point processes, in particular we give bounds on the eigenvalue repulsion, 
large gap probability, identify the limiting intensity and provide a central limit theorem. 

In the decaying model, the limiting processes are the same as the point processes obtained as the 
bulk scaling limits of the /3-ensembles of random matrix theory. In the critical model, the eigenvalue 
repulsion is much stronger. 



The critical model 

For very small values of a, this matrix behaves like the discrete Laplacian - its eigenvalues are 
locally close to periodic and its eigenvectors are extended. The discrete measure constructed by the 
square of the coordinates of the normalized eigenvector will not be concentrated on any small set 
of points. For large a, the matrix is close to diagonal, with eigenvalues dropped independently at 
random (Poisson statistics) and eigenvectors are localized. The goal of this paper is to examine the 
nature of the transition from extended to localized eigenvectors, and the corresponding eigenvalue 
statistics. 

The matrix Hn is a perturbation of the adjacency matrix of a 1-dimensional box. When the 



variance of Vk does not depend on n, eigenvectors are localized (Carmona et al, 1987 Kunz and 



Souillard 1980 Goldsheid et al, 1977) and the local statistics of eigenvalues are Poisson (Minami 



1996 Molchanov, 1981). For the perturbed adjacency matrix of higher-dimensional boxes, localiza- 



tion (Aizenman and Molchanov 1993 Frohlich and Spencer 1983) and Poisson eigenvalue statistics 



(Minami, 1996) hold if the variance is a sufficiently large constant. In dimensions three and higher. 



for a small constant variance, it is widely conjectured that one gets random-matrix type statistics 
of eigenvalues and extended eigenfunctions, while for two dimensions the opinions vary. 

Our regime, where the variance of the random variables vi are of order n^^/^ captures the 



transition between localization an delocalization. We will use the methods developed in Valko and 



Virag (2009) to analyze the asymptotic local spectral properties of if„. 



If there is no noise (i.e. a = 0) then the eigenvalues of the operator are given by 2 cos(7rA;/(ri + 1)) 
with k = 1, . . . ,n. The asymptotic density near E G (—2, 2) is given by ^ with 



p = piE) = l/y^l-EyA (3) 

which suggests that one should scale by n near E to get a meaningful limit. Thus we will study the 
spectrum A^ of the scaled operator 

pn{H^-E). (4) 

We will use the well-known transfer matrix description of the spectral problem for Hn- The 
one-dimensional eigenvalue equation Hnip = pip is written as 



where 



T{x) := I I and Mf, := T{p - Ve)T{p - v^^i) ■ ■ ■ T{p - Vi) 



Then p is an eigenvalue of Hn if and only if 

M„(j)^.(;|, (6) 

for some c G M, or, equivalently (M„)ii = 0. In view of (4) we parametrize p = E + —. We will use 
the notation M^ to emphasize dependence on A, and use the similar notation for other quantities. 
Setting 

pn y/n 
we have 

M/ = T{E + ee)T{E + e^-i) ■ ■ ■ T{E + ei) hi < i < n. (8) 

The scaling of ve = abJij ^fn ensures that, with high probability, the transfer matrices M^ are 
bounded and the eigenfunctions are delocalized. 

Theorem 1. Given E G (—2,2), R < oo and a < oo, there exists a constant c so that for all 
sufficiently large n and all t > 0, the following two statements hold with probability at least 1 — c/t. 
(1) We have 

max TiMfM}* < t. (9) 

0<^<n,|A|<i? 



(2) For each eigenvector ip of Hn, normalized by Yll=i I^^P — ^ ^^^ corresponding to an eigenvalue 



fie[E-^,E + -^], we have 



< \iJef + \iJe+if < 



2e 



0< i<n. 



(10) 



(n+l)t2 ^ '-^' ' '-^-^^' ^ n + l 

In order to understand the interaction of the eigenvalues near E and of the corresponding 
eigenvectors, one would ideally like to derive a limiting diffusion process for (jsl). The starting 
observation is that M^ cannot have a continuous limit. The obstacle is that for large n each 
transfer matrix T{E + Sk) in (tsl) is not close to / but to T{E). Thus we are led to consider, instead 



of M^, the matrices 



which will evolve regularly. 



Q^ = T-^{E)M^, 0<£<n, 




fir 




si^o loot 




Figure 1: Simulation of the entries of the first row of Qj with ii^ = 1, A = 25 and n = 10000. 



To control the correction factor T ^{E), we diagonalize T{E) = ZDZ ^ with 
D=\ ' \, Z=\ ' ' \, z = E/2 + t^l-{E/2y. 
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Theorem 2. Assume < \E\ < 2. Let B{t), B2{t) , Bs{t) be independent standard Brownian motions 
in M, VV(i) = 7f ('B2(t) + iBs{t)). Then the stochastic differential equation 



dQ' 



i\ Q \ / idB dW 

\dt+ I _ 
-i\ V dW -idB 



Z~'Q\ Q\0) = I 



(13) 



has a unique strong solution Q'^(t) : A G C,t > 0, which is analytic in A. Moreover with r = {crpy 



{Q 



A 

[nt/rj 5 



0< t< r) ^ iQ^^^it),0<t<T 



in the sense of finite dimensional distributions for A and uniformly in t. Also, for any given 
< t < T the random analytic functions Q'^nt/ri converge in distribution to Q^^'^{t) with respect to 
the local uniform topology. 

Theorem |2] is a one-dimensional version of a more general quasi-one-dimensional theorem that 



appears in Valko and Virag (2010). This proof, which predates the one in that paper, is included here 



for completeness. The preprint Valko and Virag (2010) was followed by the preprint of Bachmann 



and De Roeck (2010), who, in independent work, also study SDE limits of transfer matrices. Their 



starting point the so-called DMPK theory in the physics literature, which is essentially the study 
of diffusive limits of quasi-one-dimensional random Schrodinger operators from a slightly different 



point of view. We refer the reader to Bachmann and De Roeck (2010 ) for a discussion of this theory. 



One of the novelties of our approach is that it allows for studying the dependence on the eigenvalue 
A, which in turn allows us to deduce the scaling limit of the spectrum, the main focus here. 



The introduction of T "^{E) in Theorem ^ has the effect of changing the boundary condition for 
each n, so for the next result, we have to pass to subsequences. 

Corollary 3. Suppose that rij is a subsequence so that z^^ converges. Then T^^{E) — )■ T and 
the random matrix-valued analytic functions M^, converge in distribution to TQ'^/'^{t). Moreover, 
A„, converges in law to the counting measure of the zeros of the random analytic function A ^^ 

[TQ^/^(r)]n. 

Note that the sequence M„. has a limit if z"^^ converges. For A„, to have a limit, only the 
convergence of z"^^^ is needed. If it does, then the possible limits of T"^ are ±T and [T(5(t)]ii and 
— \TQ{t)]ii have the same zero set. 

Since if„ is symmetric, the limiting point process will live on M. The point process can be more 
effectively described by a scalar SDE. We first note that for any a,h eM? we have 



z-^ 




p I ai — bzi 
2 \ ai — bzi 



so Z ^ maps real vectors to vectors with conjugate entries. Since for A G M the transfer matrix Q^, 
and the limiting process Q^{t) will also be real, we can write 



Z-'Q\t) 



iip^{t) 
kpHt) 



for some complex numbers ip^{t) where ■?/''*' (0) = p/2 (the extra i in the above definition makes this 
and some upcoming formulas nicer). We will define the phase function (p^{t) by 



This uniquely determines ^p^{t) assuming that it is continuous in t (as long as detQ'^(t) 7^ 0, which 



follows from (13)). Ito's formula then gives an SDE for the evolution of (p (t) and we can identify 



the zeros of [TQ(r)]ii. This leads to another description of the point process limit of A^ 




/(t) 



Figure 2: The phase function ^^{t) for (t, A) G [0, 1] x [-20,20]. 



Corollary 4 (Schrodinger random analytic functions). Consider the family of SDE 's 



dip^{t) = \dt + dB + Re 



-i-p'H^dW 



^\0) = 



(15) 



coupled together for all values 0/ A G M where B and W are standard real and complex Brownian 
motions. This has a unique strong solution and for each time t the function A i— )■ y^^it) is strictly 
increasing and real-analytic with probability one. 

Moreover, for < l-El < 2 and with r = {(jp)'^ the point process A„ — arg(2;^""'"^) — tt converges 
in distribution to the point process 



Sch, := {A : <^^/"(r) G 2nZ} . 



(16) 



Remark 5. Note that the SDE (15) for a fixed A describes a Brownian motion with variance | and 
drift A. The random analytic function A i— t- v^'^(t) is given by the values of these coupled drifted 
Brownian motions at time r. In Lemma 18 we show that ip'^^'^^r) + 9 = Lp^'^^^^/'^^r) which means 
that ScK + = {X: ^^1^ ee + 2ttZ}. 

Note that the point process Sch^ is invariant under translation by integer multiples of 27r, but 
not under other translations. To fix this, we consider a translation by an independent uniform 
random variable: 

Sch; = Sch^ + f/[0,27r]. 

This version can be described through a variant of the the Brownian carousel introduced in |Valk6 



and Virag (2009) (the same is true for Sch,-, but with more complicated boundary conditions). 



The Brownian carousel. Let (V(t), t > 0) be Brownian motion on the hyperbolic plane H. Pick 
a point on the boundary dM and let x^(0) equal to this point for all A G M. Let x^{t) be the 
trajectory of this point rotated continuously around V(t) at speed A. Recall that Brownian motion 
in EI converges to a point V(oo) in the boundary dM. 

Theorem 6 (Brownian carousel description). We have 

{A : x^/^(r) = V(oo)} = Sch;. 



Section ^contains the proof of Theorem p^ and a description of the ODE for x^{t) in the Poincare 
disk model of the hyperbolic plane. Amazingly, a (less complete) connection between random 
Schrodinger operators and Brownian motion in the hyperbolic plane had been found already in 



Gertsenshtein and Vasilev (1959). 



Note that in the previous theorems we assumed < |-E| < 2. The case ii^ = is slightly different, 
but it gives similar results. Note that in that case z = i. 

Theorem 7. Let E = 0. In that case Theorem^ and Corollaries\^ and\^ hold with the following 



SDE's in place of (13) and (15): 



.A ^r^l I '^ ^ \ ^. , I *^^i ^^^2 \ \ ^_i^A 




and 



dQ=7^Z{{^ ,^ ]dt+{ ' l^ \\Z-'Q\ Q\0) = I. (17) 

-tao2 —idoi 



dif^ = Xdt + dBi + cos{(p^)dB2 + - cos{2!f^)dt, !f^{0) = 0. 



Note that since z = i we just need to fix the remainder of rij mod 4 for 2;"^ to converge and the 
parity of rij for arg(2;^"j^^) to converge. 



Properties of the limiting eigenvalue process for the critical model 



We now discuss some of the properties of the hmit process Sch,-. We describe the eigenvalue 
repulsion, we compute the intensity of the point process and then give the asymptotic probability 
of finding a large gap. We also provide a central limit theorem for the number of points in a growing 
interval. Let Sch^[a, b] denote the number of points of Sch,- in the set [a, b]. 



Theorem 8 (Eigenvalue repulsion). For /x G M and e > we have 

P{Sch,[/i, fi + e]>2}<4exp (-(log(r/£) - tY/t) 
whenever the squared expression is nonnegative. 



(19) 



Remark 9. In case of the classical random matrix models GOE, GUE, GSE the eigenvalue repulsion 

is a lot weaker: it is of the order of e'^^^ where /3 = 1, 2 and 4 in the respective cases. 



Theorem 10 (Intensity of the point process). The intensity measure A i— )■ ESchT-(y4) has density 

^p{2TTk + x), p{y) = - 

k 



Vs 






nr 



at X. This is the density of a centered normal random variable with variance |r mod 27r (a theta 
function). 

Remark 11. If the random variables ui have a a bounded probability density g{x)dx, then the 
general Wegner and Minami's estimates for random discrete Schrodinger operators ([Minami 1996 



Graf and Vaghi 2007 Belissard et al, 2007 Combes et al 2009) give 



and 



P(A„[;Ui,/i2] > 1) < Vn<y ^hWoo (/U2 - /Ui 



P(A„[^i,/i2] >2) < —na ^\\g\Lifi2 - t^i) 



(20) 



(21) 



Since we rescale the potential by y/n, both (20 ) and (21 ) diverge as n — > oo, whereas Theorems 8 and 



To] give effective bounds. Moreover, this theorem applies to singular potentials, such as Bernoulli 
random variables ±1 with probability 1/2. 

Theorem 12 (Probability of large gaps). The probability that Sch,- has a large gap is 

P(Sch.[0, A] = 0) = exp {-^(1 + 0(1)) 



where o(l) — )■ for a fixed r as A —)■ oo. 



The following theorem shows that for large A the number of points in [0, A] is close to a normal 
random variable with mean A/(27r) and a fixed variance. 

Theorem 13 (Central limit theorem). As X -^ oo we have 

(¥;0(r),y.^(r)-A)^(eo + 6,eo + 6) 

where ^0;^i;^2 o'^e independent mean zero normal random variables with variances r, r/2,r/2, re- 
spectively. 

In particular, for 6 E [0, 27r) and fc — )■ oo along the integers we have 



Sch^[0,27rA; + ^] -k 



^0+^2 + 

2tx 



^0 + 6 

27r 



(22) 



The decaying model 



The decaying model ([2]) can be thought of as the truncation of a discrete Schrodinger operator 
on the infinite half line with potential Vk = acok/vk. Similar operators have been studied in the 
literature (see Delyon et al, 1985 Kiselev et all 1998, and references therein for earlier works). In 
these works, the standard deviation of the fcth diagonal element is on the order of k^"' for a > 0. 
Depending on a, the operator has different generic spectral properties. 

• Slow decay: for < a < 1/2, the spectrum is pure point with probability one. 

• Fast decay: for a > 1/2, the spectrum is absolutely continuous with probability one. 

• Critical decay: for a = 1/2, and small enough a, the spectrum is singular continuous on an 
interval (— c, c) and pure point on (—2, — c) U (c, 2), with probability one. 

It is natural to investigate the fine eigenvalue statistics in these three cases. Motivated by this 



question, Killip and Stoiciu (2009) described the local behaviour of the spectrum in the context of 



random CMV matrices, the unitary analog of one-dimensional discrete Schrodinger operators. 
Our decaying model corresponds to the critical case. It will be convenient to reverse the indices 



to have Vk = acuk/y/n + 1 — k. We scale near E G (—2, 2) \ {0} and define Q^ as before. This 



process will converge to an SDE similar to (13), but the convergence will only hold on [0, 1) 
Theorem 14. We have the following limit on [0, 1).' 



dQ' 



1 



iX 
-iX 



dt + 



ap 



idB dW 
dW -idB 



Z-'Q\ Q\0) = I. (23) 



2 \ \ -i\ I "" ' v^r^i 

in the sense of finite dimensional distributions for X and uniformly on compacts in t. 



From (23) and Ito's formula it follows that the phase function ip'^ (which can be defined the same 



way as in (14)) will satisfy the following SDE on [0, 1): 



dip^{t) = \dt + 



ap 



dB + Re 



*'^'(*)rfw]V ip\0) = 0. 



(24) 



Note that since f^ jz~j;dt = cxd, the process f'^it) does not have a limit as t — > 1. However the 
relative phase function a^ = ip^ — (f^ will converge and its limit will describe the point process limit 
of the spectrum. 



Theorem 15. Let a^(t), <t <1, a^(0) = 0, be the solution to 



da^{t) = Xdt 



ap 



Re 



-ia^{t) _ 



l)dW 



(25) 



The function g{X) = ^^ lim(_i.i- a;'^(t) is integer valued and non- decreasing. // n — t- oo then the 
scaled eigenvalue process (see ^) converges to a point process A with counting function g. 



Applying the time change t = 1 — e ^"^^^ for the SDE (25) we get 

a^(0)=0, tG[0,oo) 



da^ = Ave"4*cit + Re 



- l){dBi + idB2) 



(26) 



where Bi,B2 are independent standard Brownian motions and (3 



(apy 



this is precisely the SDE 



that describes the SincA process given in Valko and Virag (2009). 



Corollary 16. The point process A agrees with the point process Sinews, the bulk limit of the beta 
Hermite ensembles of random matrix theory with (3 = j^^- 

The Hermite /3-ensemble is a finite ensemble with joint density 



z-,]^W\^^-H'^--^^ 



,A? 



i<i 



this suggests that the eigenvalue repulsion is of the order of e^'^^ (in the sense of Theorem^, and 



this can be proved using (25). In Valko and Virag (2009) it was proved that Sinc/j is translation 



invariant with density (27r) ^ this provides the analogue of Theorem 10 The asymptotic probability 



of large gaps was identified in Valko and Virag (2010). As A — )■ cxd we have 

'/3 1 



P (Sine/3[0, A] = 0) = (x/3 + o{\))X^' exp 



■ — A^ + 

64 



8 4 



(27) 



with 7/3 = Hf + | — 3j and < x^ < oo. 



We will also prove the analogue of Theorem 13 



10 



Theorem 17. As X ^ oo we have 

1 



^ (Sme^[0,A]- — 

A/log A V ^^ 



mo, 



/37r2 



An n — )■ oo version of this theorem for finite matrices from circular and Jacobi /3 ensembles was 
shown by |Killip| ( [2008| . 

Section [2] contains the proofs about the various properties of the limiting point processes. Section 



3 will describe some connections to the Brownian carousel introduced in Valko and Virag (2009) 



and prove the theorem about the carousel representation of the limiting point process. In Sections 
111 and Is] we will provide the proofs for the scaling limit of the spectrum for the first model (with 
the constant variance) together with the delocalization of eigenvectors. Section M will deal with the 
proof in the case of the second model (with the decaying potential). Finally, the Appendix (Section 
[7| contains the proof for the existence of unique analytic solutions for the discussed SDEs and a 
technical proposition about the convergence of discrete time Markov chains to stochastic differential 
equations. 



2 Analysis of the limiting point process 

In this section we will provide the proofs for the theorems related to various properties of our 
limiting point processes. 

We will first show a translation invariance property for the phase function ip. 



Lemma 18 (Invariance). For every 6 eM. we have 



•^ ,„A/ 



^^-^{t) + et = if\t) 



as functions of X,t. 



Proof. From (15) it is clear that v^^(t) = y^^ ^{t) + 6t satisfies the following one parameter family 
of SDEs: 



d(f^{t) = Xdt + dB + Re 



-i^'-''H)dW 



^\0) = 0. 



Since (f (t) = ip^{t) - 9t we have 



Re 



-i^'-'ii)dW 



Re 



-*^^(t)^yy; 



where W(t) = / e'^^'dW 
Jo 



is also a standard complex Brownian motion. Thus (f^ satisfies the same SDK system as ip^ with 
a different driving Brownian motion. The uniqueness of solutions shown in the Appendix implies 
that they indeed have the same distribution. D 



11 



In order to study the point process Sch^ we will use Corollary |4} Note that for A G M the 
function a^(t) = (p^{t) — '^^{t) (which we will call relative phase function) satisfies the following 
SDE system 



da^{t) = Xdt + Re 



-iaHt) _l]dZ 



a^(0) =0, 



(28) 



where Z is a. standard complex Brownian motion with dZ = e ^'^ dW. For any fixed A G M this 
can be rewritten as 

da^ = Xdt + V2sm{a^/2)dB^, a^(0) = 0, (29) 

where B^ is a standard Brownian motion with dB^ = — v2Re e~*" ^'^dZ . 

The quantity ^(^^{t) gives a good approximation for the number of points in [0, A]. Indeed, by 
Corollary |4] we have 



2tt 



a^/"(r)-Sch,[0,A] 



< 1. 



(30) 



Proof of Theorem^ (Eigenvalue repulsion). We will give two proofs of this theorem. The first one 
uses the SDE representation of Corollary |4} A second proof, at the end of Section |3} uses a geometric 
approach via the Brownian carousel. 



For yU = by ( 30 ) we have 



P(Sch^[0,£] > 2) < 



a 



'/^t) > 27r) 



(31) 



and the same holds for other yU, with a^''^ replaced by Lp^'^^^'''^ — cp^''^ , which satisfies the same SDE. 
Since this SDE is the only thing we use we can restrict our attention to yU = 0. 



Introduce the new process Y = log(tan(Q;/4)) which is well defined for a G (0, 2tt). By (29) and 
Ito's formula the process Y satisfies the SDE 



— cosh{Y)dt + - tanh(F)(it + —dB 



dY = ^^ cosh(Y)dt + - tai)h(Y)dt + ^dB (32) 

with initial condition 1^(0) = — cxd. It is clear that {a'^(r) > 27r} = {Y explodes on [0, 1]}. Consider 
now the solution Y of the SDE (32) with initial condition 1^(0) = 0. By monotonicity we have 

Y explodes on [0,r] =^ Y explodes on [0, r] =^ sup 11^(^)1 > log(r/£:). 

te[o,r] 

For e/r < 1, the inequality \y\ < log(r/£:) implies 

|£:/rcosh(y)/2 + tanh(?/)/4| < 1. 
This means that for any s > we have 



sup \Y{t)\ <\og{T/e) 
ie[o,s] 



sup \Y{t)-^B{t)\<s. 



12 



Let T be the first hitting time of log(r/£:) by Y. Then by the previous argument we have -y=\B{T) \ > 
log(r/£:) — T which leads to 

P( sup \Y{t)\ > log{r/e)) < P(^ sup \B{t)\ > log(r/£) - r) 

ie[o,T] V 2 te[o,r] 

<4exp(-(log(r/£)-r)Vr). 

Here the last inequality follows from Brownian scaling and 

P( sup \B{t)\ >x)< 2P( sup B{t) >x)= 4P(5(1) > x) < 4e~'''/l D 

te[o,i] te[o,i] 



For the proof of Theorem 10 we need the following estimate. 
Lemma 19. We have 



E {dl^\t)Y = fit) <oo. 



Proof. Differentiating (15) twice in A is justified by Theorem 24 in the Appendix. We get that for 
a fixed A, with primes denoting A-derivatives 

dip' = dt + if'dBi, dip" = ^"dBi - ^'^dB2, <^'(0) = ^"{0) = (33) 

where Bi,B2 are independent real Brownian motions with variance 1/2. This shows that the dis- 
tribution of the first and second derivatives does not depend on A, as we already know from the 
invariance Lemma [TSl 

From the first SDK we get K(p'(t) = t. Applying Ito's lemma for p''^ and (p'^ and then using 
Gronwall's inequality gives that Ey9'(t)^ and Ey9'(t)^ are bounded as functions of t only. Ito's lemma 
applied to p)"'^ gives 

W.p"{tf = [ E{p"{s)^ + p'{s)^)ds. 
Jo 

Gronwall's inequality and the fact that Kp"^{s) is bounded leads to the desired bound. D 



Proof of Theoremll^ (Intensity of the point process). We will evaluate 

r-i lim£-iE# {[<(r), ^^+"(r))] n (27rZ)} = lim(TC)-^E [Sch,[r/i, T{fi + e)]] (34) 

By Corollary the function p^{t) is analytic in A. We will use the notation w^{t) = d\p>^{t) for 



the derivative. We will first evaluate the limit in (34) by switching the interval on the right with 
[y9^(r), p^{t) + ew^ij)]. Then we will show that the error we make is asymptotically small. 
From ( 15 ) we get that xu^ satisfies the following SDK: 

dw^ = dt + Re [{-iw^)dZ^] = dt + w^lm [dZ^] , w^{0) = (35) 

13 



where dZ^ = e"**^ dW. Since the SDE in (15) has the noise term Re[(i2^'*'] and the last equation 
has Ini[(i-E'*'], the two processes are independent (for a given fixed A). From the SDEs (15) and (35) 
we get that v'^(t) is Af{fiT, |r) and Kzu'^ = r. Using the independence of (P^{t) and w^ we get 

T-^ hm£-iE# {[v9^(r), V9^(r) + ew^ir)] f] (27rZ)} = J" p{2k7r + jjt) 



e-s>0 



k& 



where p(-) is the density of A/'(0, |r). The only thing left is to show is that 



lim£-^E#{b 



^^r) + ew^ 



(r),^^+^(r))]n(27rZ)}=0. 



(36) 



We start by noting that if Z is a random variable with density /(x) and x ^^^y ^ IR+ then using 
the notation f2n{x) = Ylik(^i /(^ + 27rA;) we have 

E#{[Z + x,Z + x + i/]n(27rZ)} = / /2.(-s)rfs < |y| max /2,(s). (37) 

Jx se[0,27r) 

The same upper bound holds if y < 0. 

By (15) X = (p'^{t) + ew^{T) can be written as Xq + B{t) where ^B is a standard Brownian 
motion independent of W and Xq is measurable with respect to the a-field generated by W. Since 



the density function of B{t) mod 27r is bounded by a r dependent constant, we may use (37) after 
conditioning on W, which gives 

E#{[(p^(r) +£ro^(r),^'^+"(r))] n (27rZ)} < cE\^^'+%T) - ((/^'^(r) + £ti7^(r))| . 

To bound the right hand side, we use the integral form of the remainder in the Taylor expansion 

E|(/;^+^(r)-((/;'^(r)+£w^(r))| < ef E \dlip'' (t)] dx < ce^ 

J u 



In the last step we used the Cauchy-Schwarz inequality and Lemma 19 This proves (36) and 
completes the proof of Theorem 10 D 



Proof Theorem 12 (Probability of large gaps). Let a = a^l'' . We bound the desired probability in 
terms of phase function events: 

P (^°(r) e (0, e) mod 27r and a(r) < 27r - e) < P (Sch^[0, A] = 0) < P («(r) < 27r) (38) 

which is clear from Corollary |4] and the definition of a. To get a lower bound first note that we have 
(^°(r) = i3(T) — Re [2(t)] where Z is the driving Brownian motion in the SDE (28) for a. Since B 
is independent of Z we have 



E 



1 {v5°(r) G (0, e) mod 27r} Z >e min /2,r(a;) = ecr > 0. 

J X 
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where f2n{x) is the density of B{t) mod 27r. This means that the lower bound in (38) can be 
estimated with eCrP (o(t) < Stt — e) from below. 



Recall the SDE (29): 



da = (X/T)dt + V2sm{a/2)dB, a(0) = 0. 



In Theorem 13 of Valko and Virag (2009) the authors analyze limt_s.oo P{a^{t) < 27r) for a similar 
SDE: 

da = Xfdt + 2 sm{a/2)dB, a(0) = 

and with certain weak assumptions on / they get the asymptotics exp(— A^||/||2/8 + o(l)). The 
exact same methods with / = ^1[o,t/2] in the present case give 

P (a(r) < (2 - s)n) > exp |-^(1 + o(l))| , P («(r) < 2n) < exp |-^(1 + o(l))| • 



We omit the straightforward details. 



D 



The asymptotic gap probability for the Sine/? process (see formula (27)) was analyzed to higher 



precision in Valko and Virag (2010). Those techniques may also work here, resulting in an asymp- 



totic expansion of the gap probability. It would be interesting to see how the more precise asymp- 
totics compare to the /3-ensemble case. 



Proof of Theorem \I3\ (Central limit theorem,) 
By (15) we have ((y9''(r), ys'** 



r 



-^) = (^0 + ^1, ^0 + 6) where 



^o = B{t), Ci 



Re 



-i^°ii)dW 



6 



Re 



-*^'(*)rfW 



Clearly, ^q, C,i and ^2 are Gaussians with the appropriate means and variances and ^0 is independent 
of (^1,^2)- However, the joint distribution of (^5^2) is not Gaussian. So we need to prove that as 
A — )■ 00 the joint weak limit of (^1,^2) exists and it is given by a pair of independent normals. Let 

Z{t) = Jl^ e-'^^^'^W and Z = j^{Bi + iB2) then 



^-TfX"^- 



6 



1 r 1 r 

-. / cos{a^)dBi + -^ / sin{a^)dB2 := 6.1 + 6,2- 
/2 Jo V 2 Jo 



V2 



We will show that (^1,^2,1,^2,2) converges weakly to three independent mean zero normals with 
variances r/2,r/4,r/4. It is enough to prove that for any (01,02,0.3) G M? the random variable 
v^ = fli^i + 02^2,1 + Cf3^2,2 converges to a mean zero normal with variance r(a^/2 + 02/4 + O3/4). By 
representing the Brownian integral as a time changed Brownian motion we can see that v^ has the 
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same distribution as B (| fnio^i + 02 cos(a'*'))^ + a^ sm{a^)'^dt^ for some standard Brownian motion 
B. All we need to show is that 

2 / (oi + (32 cos(a )Y + Og sin(a; )^(it — > r(2a^ + 03 + 03) in probability. 
Jo 

Using cos(x)^ = (cos(2x) + l)/2 and sin(x)^ = (1 — cos(2x))/2 this reduces to 



cos{2a^)dt = Re / e'*" dt -^ 0, 
Jo 



cos{a^)dt = Re / e*" rft ^ 
Vo 



in probability. We work out the second claim, as the first one can be done the same way. Using 



(29) and Ito's formula we get 

1 



i\ 



die 



e'^^'dt + ^ sm{a^/2)dB^ + \ sin{a^/2fdt 

A A 



and 



e'^^dt = — (e^"^(") - 1) - — I ?,inia^l2)dB^ + - [ sm(a^/2fdt. 
iX A ./n A Jo 



'0 '^ ^ Jo 

As A — )■ 00 the first and third terms converge to a.s., while the second term converges to in L^. 
This means that their sum will converge to in probability which is what we needed to prove the 
joint limit theorem for (y9°(r), </)'*' (r) — A). 



The limit (22) follows from 
Sch,[0, 27Tk + e]-k = #{[<^°(r), <^^/"(r)] n 27rZ} - k 
where we also used that V5'^(t) is increasing in A. 



y;^/"(r) - {2k7c + e) + e 



27r 



27r 
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The proof of Theorem 17 is very similar. 



Proof 11 (Central limit theorem for Sine 13). We will consider (26) recalling that ^ is a complex 
Brownian motion with independent standard real and imaginary parts and hence 



da^ = X-e-^^dt + 2sm{a^/2)dB, a^(0) = t G [0, 00) 



(39) 



^0 



First note that a(t) = a^{T + t) with T = | log(/3A/4) satisfies the same SDE with A = 1. Therefore 

a^(oo) - a^{T) 
VMAj 
in probability. So it suffices to find the the weak limit of 

a^{T) - A 

2 TT A/log A 
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We have 



which means 



a(T) - A 



4 

1 



a{T)-X + ^^B 



2sm{a^/2)dB 



4:sm{a^/2fdt 



for a certain standard Brownian motion B. In order to prove the required hmit in distribution we 
only need to show that j^ J^ sm{a^/2)^dt — )■ | in probabihty. We have 



sm{a^/2fdt 



;log[/3A/4] 



+ 



log A Jo /3 log A /3 log A Jq 

The first term converges to 8//5. To bound the second term we compute 
4 



cos{a^)dt. 



ifiWogX 



^ /3^a-+W4 



-dt + 



Ja^+l3t/4 



log A /3A log A 



/3AlogA 



Ja^+l3t/4: 



sm{a^/2)dB 
1 



sm{a^/2ydt + 



i\ log A 



ja^+mdt. 



The integral of the left hand side is - 



il3X log A 



4ei"^mA//3 - 1 = 0((logA)-i). The integrals of the 
last two terms in the right hand side are of the order of (A log A)~^ L e^^^^dt = 0((log A)^^). Finally, 
the integral of the second term on the right has an L^ norm which is bounded by C(log A)^^. This 
means the integral of the first term on the right, (logA)"^ j^ e*" dt converges to in probability 
from which the statement of the theorem follows. D 



3 The Brownian carousel 



The SDE system (28) has a geometric interpretation using the Brownian carousel introduced in 



Valko and Virag| ( |2009| ). Recall the SDE system (|28]) 

da^{t) = Xdt + Re 



-ia^it) 



l]dZ 



a^(0) =0, 



(40) 



Here ^ is a standard complex Brownian motion. 

Consider the hyperbolic plane, let xq be a point on the boundary and let V(t), t > be hyperbolic 
Brownian motion. For a given A G M we rotate the boundary point Xq about the moving center V(t) 
with a constant angular speed A and denote its position by x'^(t). This is the Brownian carousel 



with constant speed function 1 and in Section 2 of Valko and Virag (2009) it was proved that the 



hyperbolic angle determined by the points xq, V(t),x^(t) satisfies the SDE (40). 

The evolution of x^{t) can be described by an ODE. Consider the Poincare disk model for the 
hyperbolic plane. Then the boundary points are points on the unit circle which can be described 
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by an angle x^{t) = e^^ '■*^. The hyperbolic Brownian motion V in this model satisfies the following 

SDE: 

1 - IVP 
dV = -^^dy (41) 

where 3^ is a standard complex Brownian motion. If we set Xq = 1,7'^(0) = then 

dtl' = y\ iJJ , 7^(0) = 0. (42) 

i — I V\ 

It is clear that this ODE system has a unique solution which is analytic and strictly increasing in A 
for any t > 0. Note that one usually cannot get a'^(t) from 7'^(t), however a^{t) G 27rZ if and only 
if 7^(t) G 27rZ. 

Next we will prove Theorem |6j if we add a random shift to Sch^ then the resulting point process 
can be described with a Brownian carousel. 

Proof of Theorem\^ Let U be uniform on [0, 27r] and independent of Sch^. By Corollary 111 and 
Remark |5] the point process Sch^- + U has the same distribution as the solutions of the equation 

'^^ v~) = U mod 27r. This can be rewritten as 

tt^(r) = -(/?°(r) + U mod 27r. 
Since U is independent of a'^ and ^9°, we have 

(a^(r),-¥?°(r) + [/) mod 27r = (a^(r),[/) mod 27r. 

Thus we can just look at the solutions of ^'^(t) = U mod 27r. 

For a given u G [0, 27r) the solution set of ol^{t) = u mod 27r can be described by the carousel 
construction: it is given by the set of those A G M for which the hyperbolic angle Xq, V(T), x^(T) is 
equal to u. 

By the Markov property of hyperbolic Brownian motion, the hyperbolic angle Xq, V(T), V(oo) 
is just uniform and independent of V on the time interval [0,T], so we may as well call it U . The 
claim follows. D 

We also provide an alternate proof to a version of Theorem |8] using the Brownian carousel. 

Theorem 20 (Eigenvalue repulsion). For /x G M and e > we have 

^{ScK[^,^ + e] > 2} < 4exp ( i^og{2./e)-r-lf \ ^^^^ 

whenever the squared expression is nonnegative. 



Proof. As in the first proof we can assume that /i = 0. By (31) if there are at least two points in 
[0, e] then the relative phase function must be at least 27r which means that the Brownian carousel 
had to take at least one full turn. Thus 

P {Sch^[0, £] > 2} < P {a^/^(r) > 27r} = P {7"/^ (r) > 2tt} . 



where 7 is the solution of (42). From ( |42| ) we get 

7^/^(r) < max(l - iViP)"^ = (1 - max iViP)"^ 

0<t<r 0<t<r 

which means that 

7^/^ (r) > 27r ^ 1 - — < max |Vt|^ (44) 

2n o<t<T 

In the Poincare disk model the hyperbolic distance between the origin and a point z in the unit 



disk is given by q{z) = log i j^q4 j . Thus (44) implies 



max q(Vt) > log (2Tc/e) . 

0<t<T 

The probability that the hyperbolic Brownian motion leaves a ball with a large radius r in a fixed 
time is comparable to the probability that a one-dimensional Brownian motion leaves [— r, r] in the 



same time. This follows by noting that Ito's formula with (41) gives 



dB coth(o) , 

dq = —= + y^dt 

\/2 4 

for the evolution of qiV) with a standard Brownian motion B. By increasing the drift from coth(g)/4 
to oolgg[o^i] + coth(l)/4 we see that q is stochastically dominated by 1 + tcoth(l)/4 + \B{t)\/\/2 
where B is standard Brownian motion and coth(l) < 4. Thus 

P ( max g(Vt) > log (27r/e) J < P ( max |5(t)| > log (27r/£) - 1 - r J 

< 4exp(-MV£lz^ 

which proves the theorem. D 

4 Convergence of the regularized transfer matrix evolution 

This section is devoted to the proof of Theorem [2j In order to keep the notation simple, we 
will only treat the case when m = 1, i.e. when a single value A G C is fixed. The extension to 
A = (Ai, ■ ■ ■ , Am) for TTi > 1 is straightforward. We drop A from the notation. The identity 

Tiy)T~\x) = I + { ^ ' ^ 
19 



and the recursion Mi = T(E + ei)Mi_i from dsl) implies 



= T~\E)(^^ 'Mni?)g.-i. (45) 

This shows that Qi, < i < n is a. Markov chain, the initial term Qq = I. As E is fixed, 
from now on we write T for T{E). We will work in the basis diagonalizing T i.e. we consider 
Xi = Z'^QiZ instead of Qi. (We have learned that such a change of basis has been considered 



for a slightly different problem by Schulz-Baldes (2004)). Using T = ZDZ ^ from (12), we obtain 
after simplification 

T-' f ° ^' V' = '-^ZO.Z-\ O, := ( \, '" 
\0 J 2 \^ _^2^ -1 

Therefore Xi is a Markov chain with the initial condition Xq = I and given by the recurrence 

Xi = X,_i + f/,X,„i, Ue = U^ := zp£,0,/2 (46) 



Because of the oscillating factors z'^^'^^\ the term UiXi_i is too rough to approximate a stochastic dif- 
ferential. However, on a mesoscopic scale 1 <t^ K <t^ n, the difference Xi^x — Xg = J2j=i Ue+jXi^j^i 
becomes a good approximation for a stochastic differential because the oscillations cancel in the 
sum. For convenience of the reader, we first present a heuristic derivation of the limiting SDE and 
then we give a rigorous proof. 

Heuristic proof. As Xi+K = {I + Ui+k) ••■(/ + Ui+2){I + Ui+i)Xi ~ (/ + X]j=i Ui+j)Xi, we have 

K i [ ^ \ 

Xe+K — X^ ~ 2_^ Ue+jXi = — I 2_^ ^e+jOe-i-j 1 Xg. 

We look separately at the drift and the noise contributions, i.e. we split 

j=l \ j=l / \ j=l J 

Since —2 < i? < 2, we have \z\ = 1 and z"^ ^ 1, which implies that '^j^iZ'^^'^^'^'^-'^ is bounded for 
large K. With At = K/n we then have 

_ XAt f i 
2 1^ -z 
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For the noise term we write 



u 



—lap 

2K^' 



K 







with 



K 



K 






^2l>+2j 



(47) 



3=1 

In the hmit K ^ oo, {$,i^k, ^GCe,K, ^'^Ci^k) is a mean zero Gaussian vector {$,£, Re(i, lm.Ci) whose dis- 
tribution is determined by the covariance matrix. Computing the covariance matrix is equivalent to 
computing the hmits of the expectations of ^|j^, ^e,K, Ce,K, Qk ^"^^ IO,ft:P since EReCImC = |lmEC^, 
E(ReC)^ = \ (E|CP + ReEC^), E(ImC)2 = \ (E|CP - ReEC^), E^ReC = ReE^C and E^ImC = ImE^C- 



Using (47) we get 

Ee| = E|C|| = 



K 



K 



1, niQ= hm iK-^Yz'"'^^' 



and ECl 



hm K-'yz''+^^. (4J 



The first sum in (48) converges to zero. The assumption E G (— 2,2)\{0} imphes that \z^\ = 1, 
z^ 7^ 1 and therefore the second sum in (48) converges to zero as well. Thus asymptotically C,e,K 



and (i^K are independent standard real and complex normals. Collecting our estimates we formally 
get the SDE 



dX 




Xdt + 



ap 



idB dW 
dW -idB 



X, 



X(0) 



(49) 



from which Theorem p^ would follow after rescaling time and A. In the case i? = 0, we get E(^| = — 1, 
which implies that asymptotically ^i^k and ImCe^K are independent standard normals and ReC^,^ = 
0. In this case we formally get the SDE 



dX 




Xdt + 



ap 



idBi 




X(0) 



(50) 



D 



2 \ -idB2 
where Bi,B2 are independent standard Brownian motions. 

As we will show these computations can be made rigorous. 
Proof of Theorem^ In order to make the convergence argument precise, we use Proposition 26 



which is a slight modification of Proposition 23 in Valko and Virag (2009 ). We show the convergence 



in case of a single A G C , the proof for finite dimensional marginals in A is very similar. 

We will prove that X" = Z~^Q'^Z converges to the solution of the SDE (49), from this the 



statement of the theorem follows. We can identify the 2x2 complex matrix X^ with a vector in 
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M.^ by taking the real and imaginary parts of the entries. From (46) one gets that the conditional 
distribution of X^^ — X^ given X" = a; is the same as that of 



YPix) 



iX 



iapue+i \ 



^2(^+1) 



2n 20i J \ _^2(£+i) _i 



X. 



From this b"'{t,x) and a"'(t, x) are computable. The function b"'{t,x) will be a vector in M^ corre- 
sponding to the complex matrix 



iX 



z 



2(^+1) 



.^2(^+1) 



X with 



[nt\. 



The asymptotic variance a"(t,x) is a bit more cumbersome to write down, it is an 8 x 8 matrix 
with entries which are linear combinations of terms of the form of AjAk, AjAk and AjAk with 
j,k G {1,2}^ where 



A 



lap ( 1 z2(£+i) 

_^2(^+l) _i 



with 



i = [nt\ . 



Clearly the coordinates of a"(t,x) are bilinear functions of x and x with bounded coefficients (de- 
pending on ap and various powers of z, z). 

The functions a{t,x),b{t,x) can be obtained from a'^,b"' by writing zeros in place of the (non- 



trivial) powers of z and z, these are clearly C functions. Condition (67) follows from the fact that 
n~^ sup£ J2j=i -2^"' ^^^ ^~^ sup^ X]j=i ^''"' both converge to 0. Because of this in the integrals of (67) 
the z terms will vanish in the limit and by the construction of a and b the other terms will cancel. 



Condition (68) is straightforward since b^V" are linear and a^a"^ are bilinear functions of x,x with 
bounded coefficients. The condition (69) is a consequence of the assumption E|a;^|^ < oo and since 



Xq = I the last condition is also satisfied. 



Thus we can apply Proposition 26 and the only thing left is to show that the functions a{t,x) 



b{t,x) correspond to the variance and drift functions corresponding to (49). The fact that the drift 



function agrees is straightforward. To check the variance one needs to turn (49) into a real vector 



valued SDK which basically means that we need to take independent standard real and complex 
standard normals B and W and compute the variance of the random vector corresponding to 



ap 



iB W 
W -iB 



X 



Using Ei?^ = E|iyp = 1 and "KBW = EW^ = one can check that we get exactly a(t,x) which 



finishes the proof of (49). A time-change and the reparametrization A — )■ X/t gives the SDK that 
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is independent of ap 

I ( i\ Q \ I ( idB dW \ 

and we get the claimed SDE through multiphcation on the left by the matrix Z. D 

The same argument works for the proof of the first part of Theorem [7j The only difference is 
that in that case z = i thus z'^^ = z^^ = 1 and a{t, x) will be defined accordingly. 



5 Convergence of the rescaled eigenvalue process 

In this section we prove the delocalization result (Theorem [I]) and the point process limit theorems 
(Corollary g and g. 

Tightness bounds 

Lemma 21. Let {Xk{z) : 1 < k < n,z ^ C) be random di x d2 matrices whose entries have 
finite second moments. Assume that Xk{z) is analytic in z and it is a martingale with respect to a 
filtration J^j. ■ Then for every ri < r2 < oo and t > 

-l'2~r'l i I ^^^ -^ / ^i6r2\ y ( o^^^2\* , 



¥{ max TiX J z)Xki zY >t} <t-^^—^— WIiXJe"''''^)XJe"''^''Yde 

ll<A;<n,|z|<ri J r2-ri2T[J_^ 

Proof. Since each entry Xk{z){i,j) is analytic in z, the Poisson formula and Jensen's inequality 
gives 

,2 . ^2 + n 1 r I ,^er2^'■ -M^ 



\Xkiz)it,j)\' < ^^^— / Xkie''^')it,j) de, for |z| < r,. 



Summing over all i,j gives 



TTXkiz)Xkizy < — -^^Xk, for \z\ < ri, 

r2 — Ti 2lT 



where 



Xk := r TTXk{e'''-')Xk{e''''rde 

As Xk is a submartingale, the statement follows from Doob's inequality and Fubini's theorem. D 

Proof of Theorem^ Let Sa = H = T{E+ — ). For large enough n and all complex A, |A| < r2 := 2R, 
the eigenvalues and eigenvectors of S are close to those of T{E), we can write 

Ex = ZxDxZ^' (52) 
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with \\Zx — Z\\^ , \\Dx — -DII2 , \\Z^ ^ — Z ^H^ all bounded by c/n with c depending on E, R. Since 
T[E) has unit length eigenvalues we can find another constant C = C{R, E) so that the eigenvalues 



of S^, \k\ < n are uniformly bounded by C. Using this with the decomposition (52) we get that 
TrS'^S*'^ is also uniformly bounded for \k\ < n. 



Setting Si = E Mi we have, analogously to (45), that 



^_£ —an ^I'^ijJi ^^ 
Sf. = Si^i + SiSi^i, Si := ^ \ ] ^ . (53) 

\ 

Since K£i = we get 

Taking the trace and conditioning on Si-i we get 

TiESiS; - TtESi-iS*i_^ = ETiSeSi^iS^^^S; = ETTSe-iS*i_^S;Si 

= Tt{ESi^iS;_-^){ES;Si) < TT{ESe-iS;_^)TT{ES;Se) 

In the last step we used that if A, B are positive semidefinite matrices of the same dimension then 



TyAB < T1AI1B. Using the (52) and the bounds on Z\,D\ one gets that Ti{E8l£i) < co'^jn and 
it follows that for all n > no, < £ < ra and A G C with |A| < r2 we have 

E\lx{S\Si)\ < Ci. 



Note that (S'^ : < ^ < n) is martingale analytic in the parameter A. Then Lemma 21 implies that 
(|9| holds for S*^ in instead of M/ with probability 1 — c/t. To translate the result for Mf we use 
the estimate 

TrM^M; = TrS*^S^5^5; < TrS*^S^ TrS^S; < CTiSiS;. 

To prove the second part of the theorem it is enough to show that if we assume that TrM^M^* 
is bounded by t uniformly in £ and A then ( [Io| ) holds. Let ?/' be a normalized eigenvector of Hn 
corresponding to the eigenvalue e = E + ^ E [E — -^, E + -^] and let ^i = (5,"^^)- For each 
< k,i < n, the transfer matrix description of the eigenvalue equation gives "^k = ^fc^o = 
Mfc(M^)~^\E'£, ipQ = ipn+i = 0. Since for the induced operator norm ||-||2 2 we have 

II AII2 2 = VA^ax(Av4*) < VTtAA\ 



we get the bound ||M£||2 2 < V^- But Mi is a 2 x 2 matrix and det Mi = 1 so ||M£||2 2 = \\M^ '^\ 



2,2' 



This leads to IIM„ ML^ < -\/t and 

II ^ II z,z 



2 ^ I! A# 112 II A#-l||2 lliTr 112 ^ j-2 lliTr 112 
2,2 



|^fc|l2< ||Mfc||^2P^7 22 11^^112 <i 11^ 



Summing the last inequality over all < A; < n gives 2 < (n + l)t^ jl^^jl and summing over all 
<£<n gives (n + l)||$fcf < 2^2. D 
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Proof of Corollary [3] 

By ^ and ( 11 ) for each n, the rescaled eigenvalues A^ are given by the zeros of the random analytic 
function gn '■ C — > C : 



.A/ 1 




^„(A) := det ( g;^, I ^ 1 ,5„ 1 , i?n := T-" I ^ | . (54) 

Our assumption is that along a subsequence Uj, Bn converges to a vector S G M^. It follows from 
Theorem |2] that for any fixed (Ai, ■ ■ ■ , A^) G C™, the random vector (gnji^i), ■ ■ ■ ,gnji^m)) E C^ 
converges in distribution to a random vector 



{g(\^),---,g{\J):= f det f g^^/^(r) f J j , 5 j , ■ ■ ■ , det f g^-/-(r) f ^ J,B\\ (55) 

where (g^i/^(t), ■ ■ ■ g^™/^(t)) is the solution to the SDE ([l3|. We need to show that the family of 



distributions in (55), indexed by (Ai, ■ ■ ■ , A^) G C", defines a random analytic function g{X) and 
that the mode of convergence gnjW -^ 5'(A) is strong enough to ensure convergence of zeros. 

We will use the following notions of convergence. Let A{D,'R'^) denote the space of analytic 
functions from a connected open set D in C to C'^. We equip A{D,'R'^) with the metric 

oo II „ II 

d(/,^):=^2- y ~^ll- , where \\hl := max \\hiz)\\ . 



T=\ 



Then (^A[I)^^^\d) is a complete separable metric space and convergence in d is the local uniform 
convergence. A random analytic function in ^(D,]R'^) is a measurable mapping u ^ f^ from a 
probability space {Q,J^,P) to {A{D,'R'^),B), where B is the Borel a-field generated by the metric 
d. The law of / is the induced probability measure p/ on {A{D, M.'^),Bd). A sequence f^ of random 
analytic functions is said to converge in law to a random analytic function f^ if pf^ — )■ pf in the 
usual sense of weak convergence. 

Proposition 22. Suppose 

(1) /^ is a sequence of random analytic functions in A{D, M'^) such that for every < r < oo, 

lim P \ sup |/^(A)| > C !^ = 0. (56) 

C^-oo (^£>l,|A|<r 

(2) for each m > 1 and A = (Ai, A2, ■ ■ ■ , A^) G C™" there is a probability distribution u^ on C"^ and 
the random vector (/^(Ai), fi{X2), ■ ■ ■ , fii^m)) G C™ converges in law to u^ . 

Then there is a random analytic function f^ in A{D^W'') such that f^ converges in law to f^ . 
Moreover for eachX = (Ai, A2, ■ ■ ■ , A™) G C™, (/'^(Ai), /'^(Aa), ■ ■ ■ ,/"'(Am)) G C™ /ios distribution 
u~\ 
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Proof. For each disk Dr := {A G C : |A| < r}, the bound in (56) together with Montel's and 
Prokhorov's theorems imply that a subsequence of fe restricted to Dr converges in law to a random 
analytic function fr on Dr- Then by a diagonal argument, there is a subsequence of fg^, such that 
for each integer r, the restriction of fi^ to D^. converges to to random analytic function /^ on Dr. 
The distributions of the functions fr are consistent with respect to restricting to smaller discs, and 
thus there is a random analytic function / on C such that fe,^ — ?• / in law. Condition (2) is strong 
enough to ensure that / is unique and thus fe ^ f- □ 

Let A = A{C , C ) and Aq '■= A\ {0}, i.e. we discard the identically zero function. Ai denotes 
the set of nonnegative Borel measures on C , that are finite on bounded subsets of C . We consider 
the local weak topology on A^: a sequence /i£ G A^ is said to converge to /i G C if for every 
continuous function ip : C — )• M of compact support, f ipdfii — )■ f ipdfi. For / G ^O) we denote by 
lif the zero counting measure of /, i.e. /i/ = X]f(z)=o"^('^)'^('^)' where m{z) is the multiplicity of 
the zero z. As an elementary consequence of Cauchy's integral formula, we have that for fi, f G Ao, 
d{fe, /) -^ implies /i/, -> fif. 

A random measure in A^ is a measurable function u ^ fi'^ to Ai (with the Borel a-algebra). 

If /^ is a random analytic function in A with P(/ = 0) = 0, then /xj^ is a random measure in 
A4. li fi converges in law to f^ and P(/^ = 0) = P(/ = 0) = 0, then the corresponding random 
measure fif^ converges in law to fifu^. 

We can now complete the proof of Corollary [3| The appropriate part of Theorem [7] can be 
proved the same way. 

Proof of Corollary \^ The tightness bound ([9]), Theorem |2] and Proposition 22 guarantee that the 
random analytic function C 9 A — t- Q^ G Af2(C) converges in distribution to the random analytic 
function C 9 A — )■ Q^^^('r) G M2(C) as ri — )■ oo. Then the random analytic function (?„ (A) defined 



in (54) converges to the random analytic function g(X) defined in (55). In addition it is easy to see 



that f{gn = 0) = P((7 = 0) = 0. Thus figu, converges in law to figu^. D 

The phase function 



Proof of Corollary^ The existence and uniqueness of the analytic solution of (15), as well as the 
monotonicity of v'^(t) will be shown in Section ul 

To prove the second part of the theorem we will first assume that z'^^~^^ converges to e*^. Then 
T"j (E) converges to a matrix T and 

„,, /z"^0\ ( z z\ ( ze~'^ \ f e-'^ e'^ 

TZ = limT"^ E Z = Zlim = = 

\Oz"^y \liy\0 ze'^ J \ ze-'^ ze'^ 
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By Corollary |3] we need to identify the zeros of 

[fg^(r)]n = [fZX\r)Z-%, = [fZX^r)],, 



where X satisfies the SDE (51 ). By linearity, X := XZ satisfies the same SDE, but with initial 
condition X^(0) = Z-\ 



Note that if A G M then X^ is a matrix of the form 



a b 
a b 



. Indeed, this holds for t = and 



it is preserved by the evolution by (49). So we have 



[fZX\T)l, = e-''X\T)n + e''X\Th, = 2Re [e~''X\T)^,] 



We rewrite the SDE's for the matrix entries as follows: 



2dXn = iXXiidt + iX^dS + X^dW, Xn(0) = zp/2, 
2dXi2 = i\Xi2dt + iXi2dB + Xiac/W, Xi2(0) = izp/2. 



(57) 



(58) 
(59) 



Ito's formula gives 



rfdet X^ = and 2i Im [XnXia] = det X^{t) = det Z'^ = ip/2 
which shows that X^^{t) is never equal to 0, and the phase function <f^{t) is well-defined via 

y.\0)=0. 



^i<p^{t) _ ^-^nl*^) 



^x,\ity 



Ito's formula applied to (58 ) shows that (^ satisfies ( 15 ) with —dW in place of dW. Also, the zeros of 



(57) are given by the solutions of Re[e *" *'^/^iX'^(r)ii] = 0, or, equivalently, —9 — tt/2 + (p^{t)/2 E 
ttZ. 

So far we have shown that if 2;"-^+^ — )■ e*^ then 

A„^ ^ {A : v?^/^(r) e29 + n + 2nZ} = ScK + 29 + 71 



where the last equality follows by the definition ( 16 ) and Lemma 18 It follows that 



A„ - 2 arg(2: 



n+l> 



n 



Sch^. 



Since 2 arg(2"'"''^) — arg(z^""^^) is either or 2tt and Sch^- = Sch^- + 2n the statement of the corollary 
follows. n 
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6 The limit theorem for the the decaying model 



In this section we discuss Theorems 14 and 15 The proof of Theorem 14 can be done exactly the 
same way as that of Theorem 12} Since we only need to prove the convergence in an interval [0, 1 — £] 
for a given e > 0, the fact that the coefficient of the noise term blows up at t = 1 will not cause 
any problems. 



Theorem 15 can be proved the way (26) was derived for the /3-Hermite ensemble in Valko and 



Virag (2009). The proof that we present below is not fully self-contained, we only highlight the 



main points of the arguments. 



Proof of Theorem 15_. As an analogue of the continuous time phase function we define the discrete 
phase function (f^ with the identity e*'^^' = ^^[1, 1]/X^[2, 1] and the relative phase function a^ as 
(fi^ — (fi^. Note that (p^ can be defined as a continuous function in A for any fixed n which will make 
«£ a well defined function. Equation (46) can be converted to a recursion for e*"^*^ : 



where 



%\v) 



z^'Xhz-^'v) 



%He' 



Xm 



i{l + ipet/2)+ipei/2 
1 - ipe(,/2 - ipeii/2 



(60) 
(61) 



A is an eigenvalue if 




cZ-^T-"" 




which is equivalent to e**^" = — z^"+^. The discrete version of the Sturm-Liouville theory implies 
that the number of eigenvalues in a given interval [Ai,A2] is given by the number of solutions of 
e"' = — 2;^"+^ with X G [e**^" , e**^" ] . We can also count the eigenvalues using intermediate values of 
the phase function ip^. We define p^ as a continuous function in A recursively using 



J^Q — _ 2n+2 



[in-k\ 



Then the number of eigenvalues in a given interval [Ai, A2] is given by 

#([(v^nU-^^),(^n^-^^^)]n2vrZ). 



(62) 



The main steps of the theorem are as follows. The first step is straightforward from Theorem 14 



Step 1. For every < e < 1 we have (^\n(i-£\) ^ <y^{l — s) in the sense of finite dimensional 



distributions where a^(t) is the solution of SDE (25). 

The next step shows that the relative phase function a'^ cannot change too much from n{l — e) 
to n — k. 
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Step 2. There exist a constant c > depending only on a,p and A so that for every |A| < A and 
k < en we have 



E [«(i_,j),atJ A 1] < c(E dist(af„(i_,j),27rZ) +£i/2 



n 



-1/2 



+ A;-i) 



(63) 



The proof of Step 2 can be done in a similar way as in Valko and Virag (2009). By analyzing the 



recursion (60) we can get a precise estimate on E,{af_^_^\(p^ , (p^). This can be turned into a Gronwall 



type estimate for dist(a;£,7rZ) which leads to (63). (See Sections 6.1 and 6.2 in Valko and Virag 



(2009) for details.) In order to estimate certain error terms one can take advantage of the fact that 
the rotation z'^^ in (61) has an averaging effect: 



J2^'%\ 



i=ii 



£2-1 

< C(|a<?J + ^ \ae+i - a^|) 



We would like to note that this makes our case a lot easier to deal with than the one in IValko andl 



Virag (2009) where the dependence of the oscillation on i was more complicated and needed much 



more involved estimates using harmonic analytic tools. 



The next step shows that asymptotically in the formula (62) only a^.^, 'matters'. The proof is 



analogue to the one presented in Sections 6.3 and 6.4 in Valko and Virag (2009). 



Step 3. li k = k{n) — t- 00 with k/n — )■ then (f^-k converges to a uniform random variable on 
[0, 2tt] modulo 2tt in distribution. If k is fixed then (^^ — (^^ — )■ in probability. 

Now we have all the ingredients for the proof. Suppose that we want to show that for a given 
vector (Ai, . . . , A^) we have 

the number of ev's in [0, Ai], [0, A2], . . . , [0, A^] =^ (a^i(l), . . . , a^'^(l)). 



By the previous statements we can find an appropriate sequence k = k{n) 

l,...,d) ^ («^'(1),« = l,...,ci) 



00 so that 



Ittn-fc.^ 



and (^^* — (fPf. -^ Q in probability for i = 1, . . . ,d. This means that if we apply formula (62) with 
Ai = 0, A2 = Aj then the the length of the interval will converge to a'^'(l) G 27rZ and the endpoint 
will become uniform modulo 2tt. Hence the number of lattice points in the z*'* interval will converge 
to ^Q;^'(l) which proves the theorem along the found subsequence. But the argument can be 
repeated to find a converging sub-subsequence of any subsequence, and since we always get the 
same limit this shows the weak convergence along the original (full) sequence as well. D 
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7 Appendix 

The appendix contains the proof for the existence of unique analytic solutions for the discussed 
SDEs and a technical proposition about the convergence of discrete time Markov chains to stochastic 
differential equations. 

Uniqueness and analyticity of the limiting SDE's 

Proposition 23. The stochastic differential equations (13), [T^, [P^, (figp, (23), (24) all have 



unique strong solutions which are analytic in X. Moreover the solutions of (15), (18) and (24) are 
strictly increasing in A for any positive t. 

Proof. The coefficients of these SDE's are all uniformly Lipschitz so they have unique strong solu- 
tions for any finite vector A G C*. (Note that in the decaying case one may assume t G [0, 1 — e].) 

In order to show that one can realize these solutions for all values of A together in a way that 
the dependence on A is analytic requires some extra work. 

One possibility to deal with this problem is to use the smooth dependence of the solution of an 
SDK on the initial condition. We will use the following theorem which is a slight modification of 



Theorem 40 in Protter (2005). 



Theorem 24 (Protter, Theorem 40). Let /^ : M"' — )■ M, 1 < z < d, < a < m 6e functions with 

locally Lipschitz derivatives up to order N for some < N < oo. Then there exists a solution 

X{t,u,x) to 

ft j^ pt 

i = l,...,d (64) 



xi = x,+ / f^^{x,)ds + j2 / raiXs)dB: 

Jo ^=1 Jo 



which is N times continuously differentiable in the open set {x : ({x, u) > t} where ( is the explosion 
time of the solution. Moreover the respective derivatives in x will satisfy the formal derivative of 



equation (64) 



We can encode the dependence on A in (13) into dependence on initial condition by introducing 
extra variables for ReA, ImA and the extra equations rfReA = 0, rflniA = 0. Since for any fixed A the 
SDE has globally Lipschitz coefficients we will have C, = oo. This shows that there exists a solution 



to (13) which is twice differentiable in the real variables {x^y) = (ReA, ImA). The fact that we also 



get analyticity in A G C follows from the fact that the Cauchy-Riemann equations are satisfied. 
Indeed, at time t = we have dxX{t) = idyX{t) and it can be checked that the two processes 
satisfy the same SDE which means that the previous equation is preserved. 



The same proof works for (15), (17), (18). In the case of (23) and (24) the coefficients depend on 



t as well, but introducing an extra variable for t takes care of this (note that in this case t G [0, 1)). 
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To prove that the solution if^{t) of (15) is increasing in A G M we first compute the SDE for its 
derivative. 



d{dxip^it)) =dt + Re -idxip^it)e~'^'^'UW 
For a given A the derivative solves the SDE 



dt+dx^^it)lra 



-i'/'^W^yy 



da = dt -\ — —adBx, a(0) = 
V2 



dx^^O) = 0. (65) 



(66) 



with dBx = -\/2Im 
for t > 0. (Actual 



-ivHt) dw 



and a simple coupling argument shows that this is always positive 



y, in this case one can even solve the SDE explicitly.) Similar proof works for 



(18) and (24). Note that using the carousel representation of the Section [3] one can also prove the 
monotonicity for (15) and (24). D 



Corollary 25. For a given A G M the distribution of dxf^{t) is the same as 

Jo Jo 

Proof. Using Ito's formula it is straightforward to check that the process given in the statement of 



the corollary satisfies the SDE (66). 



D 



Convergence of discrete time Markov processes to SDEs 

Proposition 26. Fix T > and for each n > 1 consider a Markov chain 

(X7gM^£ = 0... [nT\), 

with E ||X"|| < oo. For x G Mf^, let Y^{x) G M"' be distributed as the increment of X^^^ — X" given 
X^ = x. ForO<t<T andx e M^ let b'^it, x) G M"' and a'^it, x) G Mf ™(M) be defined by 

b^t, x) := nErL"„,j (x), a"(t, x) := nWY^^,,^ (x)Fl"„tj (^)^- 

We make the following assumptions. 

(1) There are C^ functions a : [0,T] x M'' ^ Mf°'(M) and b : [0,T] x M'^ ^ Md{R) such that for 

every R < oo, 



sup 

0<t<T,|a;|<_R 



{a^{s, x) — a{s, x)) ds 



+ sup 

0<t<T,|a;|<_R 



(6"'(s,x) - b{s,x)) ds 



-^0. 



(67) 



(2) For every R < oo there is a constant cr < oo such that 

||a"(t,a;) -a"(t,y)|| + ||6"(t,a;) - 6"(t,y)|| < CR\\x-y\\ , 



(68) 



31 



for all n > 1, t ^ [0, T], \\x\\ < R and \\y\\ < R. The same inequality holds for a and b. 
(3) For every R < oo there is a constant d^ < oo such that 



sup E[\\Y,^ix)f]<dnn 

0<£<n,||x||<R 



-3/2 



(69) 



(4) The initial condition Xq converges in distribution to Xq with E ||Xo|| < oo. 

Then (Xj^^pO < t < T) converges weakly in D[0,T] to the unique solution of the SDE 



dX{t) = b{t, X{t))dt + g{t, X{t))dB{t), 



X(0) = Xo, (70) 

'^ —7- Md(M.) is any C^ function 



where B{t) is the d-dimensional Brownian motion and g : [0,T] x 
with 

g{t,x)g{t,x)^ = a{t,x). 

Note: one can always take g(t,x) := (a(t,x)) but it can be useful to make other choices for 

1 /9 

which a(t, x) has sparser structure than (a(t, x)) ' and the resulting SDE has a simpler noise term. 



Proof. This is Proposition 23 in Valko and Virag (2009 ) with two small changes: there the supremum 
is for all x in (1) and (3) and the functions a, b are assumed to have bounded derivatives instead of 
being Lipschitz in x. The proof is very similar, but we include it for the sake of completeness. 

Let II ■ II oo denote supremum norm on [0,T]. For a two-parameter function / and a; G M let X 
denote the integral X f^if) = /g f{s,x) ds. We recycle this notation for a function X : [0,T] — )• M 
to write X f^xif) = /g /(s,X(s)) ds. 

Because of our assumptions on a and b the well-posedness of the martingale problem follows 



from Theorem 5.3.7 of Ethier and Kurtz (1986) (see especially the remarks following the proof), 



and even pathwise uniqueness holds. This means that (70) has a solution X with initial condition 
Xg and this solution is unique in distribution. 

Let r" = inf{t : |X"(t)| > r}. The derivation of the convergence X" ^ X is based on Theorem 



7.4.1 of Ethier and Kurtz ( |1986 ), as well as Corollary 7.4.2 and its proof. These show that if the 
limiting SDE has a unique solution (i.e. the martingale problem is well-posed as it is in our case) 



(71) 



and we have Xq =^ Xg with 








(XfenX" 


-X5,x")l(t<r;)|U 


p 

— > 





1 (Xa",X" 


-X,,x")l(t<r;)U 


p 

— > 






and 



for every £,r > sup raP(|Y£"'(a;)| > e) 

\x\<r,i 



(72) 
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then X" =^ X. The theorem there only deals with the case of time-independent coefficients, but 
adding time as an extra coordinate extends the results to the general case. 



Condititon (72) follows from the uniform third absolute moment bounds (69) and Markov's 



inequality. Thus we only need to show (71) as well as the analogous statement for a, for which the 



proof is identical. We do this by bounding the successive uniform-norm distances between 



rn,L 



^n,L 



where X"' = X^,^/^, with K = \nT/L], and X"' (t) = X'^nt\- ^^ words, we divide [0, [nT\] into 
L roughly equal intervals and then set X^' to be constant on each interval and equal to the first 
value of X" occurring there. 

If a function / takes countably many values /», then for any h we have 

i 

Since X^'^ takes at most L values, we have 

||(Xfen^Jf„,L -X6_Jfn.L)l(t < r;')||oo = | |X fen .ft^X'-^-f' 1 (^ < ^")l|oO < L SUp \\I b"-b,x\\oO = Lo{l) 

\x\<r 



by (67) where o(l) is uniform in L and refers to n — t- oo. From (68), the other terms satisfy 



(Xfe„,^,.. -Xfe.,x")l(t < r;)|U < r||(6"(-,X'^'^(-)) -&"(-,X"(-)))l(- < r, 



< c,T||X"-X"'^|U 



The same holds with b replacing 6". It now suffices to show that 



E||X"'^-X"| 



^n,L 



Esup|x;'^-x;^l </(L) 



(73) 



uniformly in n where f{L) — )■ as L — )■ oo. The left-hand side of (73) is bounded by 



Esup|x;-- J2 &n(x;^)-x;'^| + Esup|- Y, ^(^^ 

^ k=le/K\K ^ k=le/K\K 

and the second quantity is bounded by Tsup^^. \b^{x)\/L. The ffist quantity can be written as EM* 

where 

-1 



M* 



max M* 

i=0,...,L-l 



M* 



max I Mi^i I , Mi^e = Xix+e - Xix V" &" (XiK+k)- 

o,...,K-i n ^-^ 



fc=0 



Note that for each i, Mi^i is a martingale. For any martingale with Mq = we have 

I 3/2 

Emax|Mfc|^ < m^MMk - Mk-^?\Fk-^^ 

k<n 



k<n 



< cn^/^maxE[|Mfc - Mk-i\^\J^k-i]- 

k<n 
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The first step is the Burkholder-Davis-Gundy inequality (see Kallenberg (2002), Theorem 26.12) 



and the second step follows from Jensen's inequality. Therefore (69) implies 

E[|M;| Vii] < c{n/Lf'^n^ = cL^ , 
which gives the desired conclusion 



L-l 



{EM*f < E{M*f < E^{M*f < cL^. 



i=0 



Letting first n — )■ oo and then L — )• oo gives (73) and (71). 
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